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Abstract 



Recently a new dark matter candidate has been proposed as a consequence of 
universal compact extra dimensions. It was found that to account for cosmological 
observations, the masses of the first Kaluza-Klein modes (and thus the approximate 
size of the extra dimension) should be in the range 600-1200 GeV when the lightest 
Kaluza-Klein particle (LKP) corresponds to the hypercharge boson and in the range 
1-1.8 TeV when it corresponds to a neutrino. In this article, we compute the elastic 
scattering cross sections between Kaluza-Klein dark matter and nuclei both when 
the lightest Kaluza-Klein particle is a KK mode of a weak gauge boson, and when it 
is a neutrino. We include nuclear form factor effects which are important to take into 
account due to the large LKP masses favored by estimates of the relic density. We 
present both differential and integrated rates for present and proposed Germanium, 
Nal and Xenon detectors. Observable rates at current detectors are typically less 
than one event per year, but the next generation of detectors can probe a significant 
fraction of the relevant parameter space. 



1 Introduction 



Recently a new type of dark matter candidate has been proposed, consisting of stable 
Kaluza-Klein (KK) modes of ordinary standard model particles allowed to propagate 
in one or more compact extra dimensions 0]. Under reasonably mild assumptions 
about the nature of the UV completion, a five dimensional (5d) theory with orbifold 
boundary conditions has a lightest Kaluza-Klein particle (LKP) which is stable, and thus 
would be present as a cold relic in the universe today. If one further assumes that the 
LKP is a neutral, non-baryonic particle, it has all of the properties required of a well- 
motivated weakly interacting massive particle (WIMP), and what remains is to determine 
the relevant range of masses and other parameters needed in order to correctly account 
for cosmological observations, and to determine the sensitivity of current and future dark 
matter searches to detect it either directly or indirectly. 

Models with compact extra dimensions possess a very rich phenomenology and have 
attracted much attention lately. Fields which are permitted to propagate in the extra 
dimensions appear to a low energy observer as a tower of increasingly massive particles, 
each with identical spin and charge. The massive states are in fact modes of the fields 
which carry (quantized) momentum in the extra dimensions, and thus the spacing of the 
tower is roughly the inverse size of the extra dimensions. The ordinary particles 
of everyday experience are identified with the "zero modes", carrying zero momentum 
in the compact directions. Models with extra dimensions may be classified according to 
which fields are allowed to propagate in extra dimensions. In brane world models, only 
gravitational fields can propagate in extra dimensions and therefore have KK excitations. 
There are other interesting models where standard model gauge bosons live in the bulk as 
well while standard model matter fermions are still confined on a 3-dimensional surface. 
And finally, there exist models -the so-called models with Universal Extra Dimensions 
(UED) IQ- in which all standard model fields have KK modes. These three classes of 
models have very distinct phenomenology. UED models are the only models with extra 
dimensions to have a stable KK particle because of a discrete symmetry (KK parity) 
which is remnant of the higher dimensional Poincare invariance. Under KK parity, all 
even mode-number particles (including the zero modes) are even, while the odd modes 
are odd. This has the desired consequence that the lightest odd mode must be stable. If 
this lightest KK Particle is neutral and weakly interacting, it provides a cold dark matter 
candidate. 

Having fermions in extra dimensions requires further assumptions about the nature of 
the compact dimensions. Generally, from the four dimensional point of view, a fermion 
in higher dimensions will not have chiral interactions, a necessary ingredient of any elec- 
troweak theory. This may be addressed in the UED context by imposing orbifold boundary 
conditions which project out the zero modes of the unwanted degrees of freedon responsible 
for vector-like interactions. For five dimensions, this can be accomplished by identifying 
the coordinate of the extra dimension y with —y, folding the circle onto the line jZ^- 
In six dimensions, starting from a 2-torus T^, one can identify either points related by 
a rotation of 180 degrees {T'^ / Z2) or by 90 degrees (T'^/Z^). One can similarly consider 
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cases in even more dimensions, but we will restrict our discussion, for simplicity, to the 
five or six dimensional cases. Having imposed the orbifold in order to recover a chiral low 
energy theory, it can be shown that the needed boundary conditions imply that there 
are terms in the Lagrangian which live on the fixed points of the orbifold transformation. 
In five dimensions, these are the points on the boundaries of the extra dimension. These 
boundary terms cannot be computed in terms of other parameters without knowing the 
UV completion of the theory. Consequently they must instead be treated as parameters 
of the UED model. We expect that the masses of the first level KK modes should be of 
order 1/-R, but will have corrections from the boundary terms which will in general be 
different for different fields. 

The most interesting cases for dark matter are when the boundary terms are such 
that the LKP is a KK mode of either a neutrino, a neutral Higgs or a neutral weak gauge 
boson. In this work we focus on the neutrino and gauge boson possibilities. Because 
of electroweak symmetry-breaking, the KK towers of the hypercharge boson B and the 
neutral SU{2) boson mix. The mass matrix for the first level KK modes (in five 
dimensions) may be expressed in the {B^^\ W^s^^) basis, 

where R is the size of the extra dimension, v is the Higgs vacuum expectation value, gi 
and (72 are the gauge couplings, and 5M^ are the boundary terms. If the boundary terms 
are induced radiatively, they should be proportional to the gauge couplings and l/E?. 
Thus, for 1/R ^ v, the matrix is rather close to diagonal, and since gi < g2 we can 
expect that the lighter particle is well-approximated as being entirely B^^\ Thus, in this 
case the LKP is a massive neutral vector particle which couples to matter proportionally 
to gi times the hypercharge. 

In Ref. m, we determined the relic density for the LKP when it is either a KK mode 
of a neutrino {i'^^^) or of a neutral gauge boson (B^^^). As seen above, these represent 
natural candidates when terms confined to the orbifold fixed points are taken to arise 
radiativelyQ as opposed to being present at tree- level 0. A variety of co- annihilation 
channels were included, with a range of mass splittings between the LKP and heavier first 
tier KK modes, and the conclusion is that in order to correctly account for the observed 
density of dark matter, the LKP masses should lie in the ranges 600 to 1200 GeV for 
5(1) and 1000 and 1800 GeV for We also noted that for a six dimensional (6D) 

orbifold T^/Z2, these mass ranges are lowered by approximately a factor of ^/2- Under 
our assumption of small boundary terms, the LKP mass corresponds to the inverse radius 
of the compact dimension and we expect all first level KK modes to have masses of this 
order. The range relevant for dark matter is particularly tantalizing because it lies just 
above the current bounds from high energy colliders Given the very large number of 

^B'^^^ is indeed the LKP when assuming vanishing boundary terms at the cut-off scale. Note that this 
kind of prescription is similar to the choice of universal soft SUSY breaking masses at the GUT scale in 
susy models. 

^ The possibility of an extra dimension at a TeV was first examined in |^ . 
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currently running or planned experiments devoted to both direct and indirect searches for 
WIMPs, the detectability of KK dark matter is an interesting question. Indirect detection 
issues have recently started to be investigated P,|lO|,p!l|. 

Direct detection of a weakly interacting massive particle typically involves searching 
for the rare scattering of the WIMP with a nucleus in a detector. As a result of the 
interaction, the nucleus recoils with some energy, which can be read out as a signal . 



The distribution of recoil energies is a function of the masses of the WIMP and the nucleus, 
and (because the scattering length for heavy WIMPs is typically of the same order as the 
size of the nucleus) the nuclear wave function. The lightest supersymmetric particle (LSP) 
is a typical Majorana fermion WIMP with mass on the order of 100 GeV, and theoretical 
predictions for its interactions at modern dark matter detectors have reached a high level 
of sophistication. In Ref. |]^, computations for the cross sections of B^^^ scattering with 
nucleons were performed and prospects for direct detection at present experiments were 
presented. In contrast with Supersymmetric WIMPs, predictions depend only on three 
parameters: The LKP mass, the mass difference between the LKP and the KK quarks 
(assuming all flavors and chiralities of first level KK quarks are degenerate in mass) and 
the "zero mode" Higgs mass. It is clear that present experiments can only probe KK 
masses below 400 GeV as soon as the mass splitting between the LKP and KK quarks 
is larger than five percents (as is found in |^). On the other hand, masses below 300 
GeV are already excluded by collider constraints And in any case, masses below 400 
GeV are in conflict with the mass range predicted from our relic density calculation |]T|]. 
Therefore, we wish to investigate detection prospects for masses above 400 GeV and ask: 
Which planned experiment will be able to probe the relevant parameter space of the 
LKP? To answer this question, we need to go beyond the wimp-nucleon cross section 
calculation and compute the event rate for a given detector. This requires to include 
the nuclear form factor. In the current article we expand upon the results of Ref. |0, 
deriving realistic estimates for event rates at modern dark matter detectors, including 
nuclear wave function effects and examining differential rates in the nuclear recoil energy 
as well as integrated ones. We find that the event rate is somewhat smaller than for the 
usual LSP neutralino WIMP with mass around 100 GeV and that in order to see at least 
several events per year, heavy (> 100 kg) detectors are needed. 

This article is organized as follows. In section II we review the kinematics of direct 
detection of WIMPs. In section III, we present the analysis in the case where the LKP 
is the first KK state of the neutrino, finding that it should most likely have already 
been observed by CDMS or EDELWEISS, and thus is excluded. Section IV is devoted 
to the more interesting case of B^^\ Our predictions for the differential and integrated 
event rates expected in Germanium, Sodium-Iodide and Xenon detectors are presented 
in section V. We reserve section VI for our conclusions and outlook. 
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2 Kinematics of WIMP Detection 



In this section we briefly review the general kinematics of WIMP-nucleus scattering. The 
number of events per unit time and per unit detector mass is, 

where m is the WIMP mass and p its mass density in our solar system^, and M is the mass 
of the target nucleus. f{v) is the distribution of WIMP velocities relative to the detector, 
fi = mM/ (m + M) is the reduced mass, q'^ is the momentum transfer four- vector whose 
magnitude is |qp = 2/i^f^(l — cos^) in terms of 6, the scattering angle in the center of 
momentum frame. |qp is related to the recoil kinetic energy Er deposited in the detector 
(in the lab frame) by Er = |qp/2M. For m M as is the case for LKP WIMPs with 
masses of order 1 TeV, Ej, is typically 30-50 keV depending on the nucleus target (but 
it can be much larger for WIMPs with velocities close to the galactic escape velocity). 
Eq. (^ may be thus rewritten 

dR 2p da 
dEr = m 

in which |qp should be regarded as a function of Er as indicated above. The differential 
cross section can be expressed in terms of the cross section at zero momentum transfer 
(To times a nuclear form factor p^ , 

d\q\^ Ap^v^ ^ ' 

where F^(|q|) is a function normalized to one at |qp = which includes all relevant nuclear 
effects and must be determined either directly from measurements of nuclear properties 
or estimated from a nuclear model, and o"o contains the model-dependent factors for a 
specific WIMP. The rate is obtained by integrating over all possible incoming velocities 
of the WIMP: 

dEr 2mp? ^' '^i^>„,„ V ' ^ ' 

where v^ax — 650 km/s, the galactic escape velocity. To determine f^m we use the 
relation between the WIMP energy E and the recoil energy Er 

1-cos^ 4mM 4mM ' ^ ' 

Vmin = {2Emin/mf^'^ = C^) 



■^In our numerical results, we identify p with the local (in our galaxy) dark matter density of canonical 
value p ^ 0.3 GeV/cm"^ ||l3|] i.e. we assume a homogeneous halo of our galaxy. This assumption is not 
obvious. In addition, estimates of the local density of dark matter in our galaxy is subject to considerable 
uncertainty and model-dependance. We should then keep in mind that local over- or underdensities could 
easily change the expected detection rate by a significant amount. 
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Assuming a Maxwellian velocity distribution for the WIMPs and including the motion of 



the Sun and the Earth one obtains [ p!^ , 



fiv] 



dv 



2VE 



erf 



+ Ve 



erf 



vo 



where ve is the relative motion of the observer on the Earth to the sun (and thus shows 
an annual modulation), and vo is the mean relative velocity of the sun relative to the 
galactic center^ Thus, the final formula for the measured differential event rate is. 



dR 



dEr AvEmfi' 



Ve 

erf — erf 

\ Vq ) \ Vq 



(9) 



The total event rates per unit detector mass and per unit time will depend on the range 
of energies to which the detector is sensitive. Thus, the actual observed rate, modulo 
experimental efficiencies, will be given by dR/dE^ integrated over the appropriate range 
of energy for a given experiment. 

Our task will now be to compute (Jq and to combine it with the correct form factor 
F^(|q|) in cases where the WIMP is z/^^^ or B^^\ To compute ctq, we must evaluate the 
effective WIMP interaction with nuclei by evaluating the matrix elements of the nucleon 
operators in a nuclear state. This in turn is determined from WIMP interactions with 
quarks and gluons evaluated in nucleon states. Traditionally, one differentiates between 
two very different types of WIMP-nucleon interactions- spin- dependent interactions and 
scalar interactions. 

Scalar interactions are coherent between nucleons in the nucleus, and the form factor 
is thus the Fourier transform of the nucleon density. The commonly used form (iden- 
tical, in the limit of low momentum transfer, to the one derived from a Woods-Saxon 
parametrization of the nuclear density |I^pl 



IS 



F\H) 



(10) 



where Ri = ^ R? — and i? ~ 1.2 fm A^/^ with A the nuclear mass number, s ~ 1 fm 
and ji is a spherical Bessel function. 



sin[grn] - qr^ cos\qr^ 



An axial-vector interaction leads to interactions between the WIMP spin and nucleon 
spin. In this case one must evaluate the matrix elements of nucleon spin operators in the 
nuclear state. The form factor is typically written as |T3,|T3| , 



FMq\) 



sio) ■ 



(12) 



^Together with p, note that vq is another crucial quantity in both direct and indirect methods of dark 
matter detection which is subject to significant uncertainties. 
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where, 



^(|q|) = alSoo{h\) + alS^M) + aoa,Sol{\q\), (13) 
ao = ap + a„, ai = Op - a„, (14) 

where the first term is the iso-scalar contribution, the second one is the iso-vector contri- 
bution and the last one is the interference term between the two. The Sij are obtained 
from nuclear calculations. Op and a„ reflect the spin-dependent WIMP interactions and 
average spins for neutrons and protons in the nucleus and will be defined below. 

3 Direct Detection of u^^^ 




Figure 1: Leading Feynman graph for effective z/'^^^-quark scattering through the exchange 
of a zero-mode Z gauge boson. 

For our first example, we consider the KK neutrino, This is almost a case which 
has been considered previously [0, the only difference being that the KK neutrino has 



vector-like weak interactions. In the non relativistic limit where <^ m|, we have an 
effective four-fermion contact interaction (Fig. |I]), 



Kl'^u^] {{Qr + 9l) [(l{x)l,,q{x)] + {g]i - gl) [q{x)i^im{x)]) , (15) 



4 sin^ Ow^m'^r 

where we have explicitly included the Z couplings to 

9r =9l =Tr—n 

2 sm 6w cos Ow 

and the gl^ and g'^^ are the left- and right-handed quark interactions with the Z boson, 

5^ = Tl-g^sin^^^. (17) 

Thus we see that the effective interaction includes both a coupling to the vector and the 
axial vector quark currents. When evaluating the WIMP-nucleon cross section, this will 
be summed over all flavors of quarks and will involve matrix elements (qj^q) and (jl'^f/y^q) , 
where the expectation values are to be understood as refering to nucleon states. 
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The WIMPs are highly non-relativistic, and thus only the time-component of the vec- 
tor u^'j^Uu is appreciable. However, the expectation value (^7075^') — |jl5|], and we are 
left with only the time-component of the vector interaction. This illustrates the predom- 
inant difference between u^^^ and a typical massive Dirac neutrino WIMP - the absence 
of spin-dependent interactions. However, since the spin-dependent contribution is usually 
sub-dominant to the scalar interaction, the resulting cross sections remain comparable. 

At the quark level, the effective interaction has the form. 



where, 



bq [uu^'^Uu] [g(a;)7^g(x)] 



'1^ 



4 sin Ow'mw 



Tl - 2Qq sin^ 9 



w 



(19) 



The matrix element (qj^q) = (q^q) simply counts valence quarks in the nucleon, and so 
the nucleon WIMP couplings are. 



G 

bp = 2bu + bd = ^ (1 - 4sin2 Ow 



V2 



2bd + bu = - 



G p 

71' 



w 



(20) 
(21) 

1/4 



for the proton and neutron, respectively. The numerical accident that sin^ 6 
renders the coupling to protons very small. The vector interactions are coherent, and 
thus we have for the WIMP-nucleus coupling, bjq = Zbp + {A — Z)bn. Thus, 



fi^G] 
2n 



:i-4sin2 9w)Z-{A-Z) 



[22] 



and the form factor entering in the differential cross section da/dc^ is given by ([T0|) . 

It is well-known that the mass of Dirac neutrinos is strongly constrained by elastic 
scattering experiments such as CDMS [16] and EDELWEISS |17]. The exclusion plots 
are presented in the m-a^ plane where m is the mass of the dark matter candidate and 
(7„ is the scattering cross section per nucleon. It is related to ctq by. 



mz 



(23) 



being the mass of the nucleon. For Ge, and m ^ M, we find ~ 2 x 10 



-39 



cm2 ~ 2 X 10-3 pb. Given that CDMS and EDELWEISS did not see any events, a WIMP 
with this cross section must have a mass > 50 TeV. This means that in order to have 
escaped detection, u^^^ would have to have masses more than ten times larger than the 
range of masses for which result in the correct dark matter relic density. While one might 
imagine that coannihilation in various channels could push up the favored u^^^ masses by 
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a few TeV , it seems unlikely that the relic density calculation could favor masses above 
10 TeV. 

To conclude this section, the KK neutrino seems to be ruled out as a dark matter 
candidate at least in the minimal UED model where the mass window prediction from the 
relic density calculation is in conflict with direct detection experiments. Let us therefore 
now concentrate on the S(^) LKP candidate. 

4 Direct Detection of B^^^ 

I 

I 

q ^ q 



Figure 2: Leading Feynman graph for effective i?*^^^-quark scattering through the exchange 
of a zero-mode Higgs boson. 



B'^^^ can interact elastically with a quark by exchanging a KK quark in the s- and 
t-channel or by t-channel Higgs exchange. The amplitude for scattering between quarks 
and B^^^ mediated by Higgs exchange (Figure ||) is, 



M, 



2 ml 



(24) 



where are the B^^^ polarization vectors, q{x) is a quark field and there are separate 
couplings 7g for each flavor of quark, gi is the hypercharge coupling, and Yh = 1/2 has 
been explicitly included in the result. We have taken the non-relativistic (NR) limit for 
the WIMPs in which we are justified in dropping tiny terms of order {ps — PbY /nT-h- The 
factor of niq in 7^ is a direct consequence of the fact that zero mode quark masses result 
from the quark couplings to Higgs, after electroweak symmetry breaking. 

We now consider the KK quark exchange, with Feynman diagrams shown in Figure |^. 
Recall the coupling B^^^-q^^^^j^yq involves a right (left)-handed projector and both g^-* and 

q]/^^ can be exchanged and will typically have somewhat different masses. Thus each 
Feynman graph of Figure ^ is actually two separate graphs with g^^'* and gj^'* exchanged. 
The amplitudes corresponding to the two diagrams of Figure 01 are: 



R/L 



q{x)YPR/L 



{i>q- i>B> + mq(l)) 



,(1) 
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,(1) 




Figure 3: Leading Feynman graphs for effective i?*^^)-quark scattering through the ex- 
change of a KK quark. Both q}^ and should be understood in each graph. 



R/L) 



q{x)YPR/Lj ■ ^ -2 — 7 PR/Lq{x) 



) 



e;(Pi?')e.(Pij),(26) 



where Yr/l are the hypercharges for the right- and left-chiral quark q. In the nonrela- 

^ [mi 



tivistic hmit pb ~ Pb' ~ {^bw^ O5 O5 0) so that M.^/^ = A^f^^ + can be rewritten: 

(27) 

Pr/l q{x) , 



_^R/L 



q[x) 



Qr/l fJ-^ 



[rriBW + Eg 



m 



(1) 



(1) 



where we have neglected the 3-momentum of the quarks in the nucleon and written 
Pq ~ {Eg, 0,0,0). Eg < 1 GeV ^ is the energy of a bound quark in the nucleon. 

We now expand this expression up to linear order in x = Egm^n) / {m?^^^) —rn?^^) ) to obtain 



-Mf ^ = -^^l{PB') eAPBM 



X] 



SgE^"" + AgE^"" 



Pr/l q{x) 



{2i 



where the coefficients Sg (scalar contribution) and Ag (spin-dependent contribution) are 
defined in equations (p^ and 

(29) 



E^""" = Yl^Y + 7-^707'^ , E^'^ = Yl'^Y - Y-f^Y = 2ze°'^'^''7p75. 



In the non relativistic limit E^^" leads to scalar interactions whereas E'^" leads to spin- 
dependent interactions. 

We will assume that all flavors and chiralities of flrst level KK quarks are equal and 
parameterize their masses by A = {nigii) — uibw) /ttibw- Summing ^Ah and J^g + J^g 
we obtain: 

(30) 
(31) 

(32) 



{M) = -it^ipB) e,.ipB') [(7, + Sg)g^''{qq) + AgiqE'^'^q) 
I {Mg)f = 4(7, + Sgfmi' + 2Al I (g7Sg) 



with, 



We sum over the different quark contributions to obtain the matrix element in a nucleon 
state. At this stage, as a first-order evaluation, we will make the assumption Eg ^ rriq. 
We recognize that the assumption that the light quarks in the nucleon are on-shell is 
questionable and that a more accurate treatment would be desirable[|. 

For the spin matrix element only the light quarks u, d, s contribute while for the scalar 
matrix elements there are also contributions from heavy quarks c, b, t 



=-p,nE^/r (33) 

Because of this distinction, we drop terms in Aq which are proportional to any power of 
the zero mode quark mass since these are negligible. Note that under our assumption both 
7g and Sq are proportional to the quark mass m^. Thus, given the normalization of the 
matrix elements /y^, each flavor contributes to the scalar interaction proportionally to its 
contribution to the nucleon mass. For heavy quarks q = c,b, t, the contribution should in 
fact be considered to be induced by the gluon content of the nucleon, with the heavy quark 
legs closed to form a loop. In the Higgs exchange case, the mapping from the tree graph 
with heavy quark external legs to the loop graph with external gluons is straight-forwardly 



handled by the formalism of [0. For the KK quark graph, as emphasized in [|r^, this 
mapping is generally unreliable because of the presence of the heavy KK quark in the loop 
with mass ~ rn^(i). Thus, we include q = c,b,t in ■jq, but to be conservative not in Sq. 
From a practical point of view, the loop-suppression renders the contribution from the 
heavy quarks irrelevant compared to the strange quark contribution, so the final results 
are insensitive to this choice of procedure. 

The coefficients Aq and Sq may be extracted from the matrix elements, 

A = '-""^T^, (34) 



ba — -r^al n n ^(mr,(i)+m^ 



1 ~ ~^1~2 _^2 \2y"BW^'"qW, 



The total amplitude squared in a nucleus state reads: 



\{M)\' = \{M)\Liar + \{M)L^, (35) 



where, 



\{M)\Liar = ^rnUZff^ + {A-Z)ff'rxFl{\ci\), (36) 
K-^)tn = f^ W^r^'""!? .. AV(J + l)xF.^^(|q|), (37) 



^ Note that this problem did not arise in the clastic scattering of a neutralino because neutralinos 
can only exchange bosons when interacting with spin-1/2 quarks. On the other hand, our bosonic dark 
matter candidate can exchange fermions with the quarks. This interaction involves fermion propagator 
and therefore involves the energy, and not only the mass, of the quark. Possibly the use of a constituent 
mass might be more appropriate in this case. 
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so that the corresponding cross sections at zero momentum transfer, ao are, 

2 

scalar '^N ( ry rB^ , (a ry\_eBW\2 



A-rr(^_,., _L ^.,^2 



47r(m5(i) + uin] 



{zfr +{A-z)fr)\ (38) 



where, 
and. 



spin 2 2 4 ^J{J + 1) /„„^ 



^ _ 0'p{Sp) + CLn{Sn) ^^^^ 
«p(n) = E + (41) 



Using A^^^^M = A(p)rf = Arf, A(")c/ = A^^^n = Am and A^*").? = A^f^s = As, we have, 

ap = i^A^i + ^(Arf + As), (42) 

an = ^Ad+^{Au + As), (43) 
-3d 6b 



where Am = 0.78 ± 0.02, Ad = -0.48 ± 0.02, and As = -0.15 ± 0.02 0. For the spin 
form factors we will also need the iso-singlet and iso-vector combinations, 

ao = ^(An + Arf) + ^As, (44) 

ai = -(Am -Ad). (45) 
3 

For Germanium and Iodine nuclei we have, assuming the WIMP mass is 1 TeV, A = 
(mg(i) — m^ii)) /m^ii) = 15%, and the Higgs mass is rrih = 120 GeV, 

Germanium : a^^"'''^ = 1.6 x 10"^ pb, ctq^^" = 1.4 x 10"^ pb, (46) 
Iodine : a^""'"" = 1.34 x 10"' pb, a^" = 1-8 x lO'^pb. (47) 

These values are somewhat smaller than what one would typically expect for neutralino- 
nucleus elastic scattering. In that case, one finds 



^scalar ^ J^^Z + {A - Z) f^f , (4J 

32 



TT 

SO that. 



~ -GVaV(J + 1), (49) 



^scalar 

^o,bw 




1 


(fr 


^scalar 


16m^(i) 
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(50) 



Using fp = mpY.fTjq/mg where f^'^^ ~ 7^ and ~ 92Thoo hhqq/2ml (note that 
and have different dimensions) where T/ioo and /i/^gg are Higgs-neutrahno-neutrahno 
and Higgs-quark-quark Yukawa couphngs (which can be found, for instance, in Ref. [|14|), 
we have, 



a, 



0", 



scalar 
0,BW 
scalar 
0,X 



91 



mw 



10" 



(51) 



Therefore we expect a^^^''"'^ to be smaller than cTq^^'ij', however the ratio generally depends 
on the precise neutralino couplings, which are complicated functions of SUSY parameter 
space. We now compare spin- dependent cross sections: 

' 1 / \4 4 

I f n-, \ m,, . 

(52) 



spin 



9i o,p \ 1 9i\ 
: — cx; — I — — I — = ~ 



spm 



48 



a$ j G^m^d) - m (i))2 



•^q^^ij \92/ V-bW -qy 

We again have a large suppression factor due to the large WIMP mass unless rrigd) is 
nearly degenerate with rriBm ■ 
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Figure 4: Spin-dependent and spin-independent WIMP-nucleon cross sections as a func- 
tion of the WIMP mass for (top to bottom) A = (mq(i) — m^(i))/m^(i) = 5, 10, 15% and 
nih = 120 GeV. 



It is common for dark matter search experiments to express their constraints in terms 
of effective WIMP-nucleon cross sections. 
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(53) 
(54) 



For 1 TeV WIMP mass, typical values are a^™^"'' ~ 10"^° pb and a^^j" ~ 10"^ pb. (For 
comparison, nucleon-neutralino cross sections are in the range 10~^^ — 10^^ pb for scalar 
interactions and 10~^ — 10~^ pb for spin-dependent interactions). ^From Fig. ^ we see 
that the cross sections may vary upward by about one order of magnitude if m^{i) is at 
the lower end of its favored range, 600 GeV, and by two orders of magnitude if in addition 
and g'-^-* are more degenerate, A ~ 5%. The dependence on the zero-mode Higgs 
mass is presented in Fig. |^. Note that theories in which the top and/or bottom quarks 
propagate in extra dimensions I^T] generically have additional contributions to electroweak 
observables through the oblique parameters S and T [|], and thus the preference in the 
precision electroweak data for a light SM-like Higgs may be misleading in theories with 
universal extra dimensions. Thus, we consider a wider range of Higgs masses than one 
would naively expect from the electroweak fits. Finally, in Fig. ^, we show a scatter plot 
of spin-dependent and spin-independent cross sections, varying 600 GeV < m^w < 1200 
GeV, 5% < A < 15%, and 100 GeV < < 200 GeV. 

In any case, these cross sections are below the reach of any currently running experi- 
ment. However, larger mass detectors composed of heavier nuclei and improved efficiencies 
will most likely change this situation in the foreseeable future. Since precise event rates 
will depend on experimental issues such as efficiencies and background rates and rejection, 
it is important to include nuclear effects in the theoretical predictions, and worthwhile to 
study kinematic distributions such as dR/dE^.. 
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Figure 5: Scalar WIMP-nucleon cross section as a function of the Higgs mass for m^(i) = 1 
TeV and (top to bottom) A = (mg{i) — m^{i))/m^(i) = 5, 10, 15%. 
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Figure 6: Predictions for iJ'^^^-nucleon cross sections in the spin-independent - spin- 
dependent plane. We have varied three parameters: m^(i) in the 600-1200 GeV range, A 
in the 5-15 % range and rrih in the 100-200 GeV range. 
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5 Differential and Integrated Event Rates 

From Eq. the number of events per kilogram of detector per keV per day is propor- 
tional to, 



dR o-Q 2 

OC ^ 

dEr mfi'^ 



« ^i^'(lql)- (55) 



Larger m and smaller ctq combine with a suppression from F^(|q|), making the event rate 
quite low. The rates are further suppressed by nuclear form factors which drop quickly 
as the recoil energy increases. 

Nucleus-A Tj^ical Recoil Momentum Typical Recoil Energy 

q = fiv |q|72M 

Silicon-23 22 MeV 11 keV 

Sodium-29 28 MeV 14 keV 

Germanium-73 68 MeV 32 keV 

Iodine-127 112 MeV 50 keV 

Xenon-131 116 MeV 51 keV 



Table 1: Typical nucleus recoil momenta and energies after scattering with a WIMP with 
mass m ^ M. 

In Table ^ we list some typical recoil momenta and energies (corresponding to a WIMP 
mass of 1 TeV and velocity of f ~ 220 km/s ~ 10~^c) scattering from various nuclei. Note 
that there is effectively a maximal recoil energy which is roughly 16 times the typical recoil 
energy listed in the table, because the maximum velocity is approximately the galactic 
escape velocity, Vesc = 650 ±200 km/s ~ 2v and qmax = '^I^Vmax- However, at such energies 
the nuclear form factor itself already provides a high suppression in the differential rate, 
such that one arrives at a good approximation to the integrated rate by integrating up to 
an energy which is four times the typical energy. Thus, it is enough for our purposes to 
present dR/dEr over a 200 keV range of recoil energy. Experimentally, it may be useful to 
look at energies > 100 keV for which we expect the background to fall off. We illustrate 
the importance of nuclear effects in Fig. |^ where we plot the scalar and spin form factors 
for il27, as a function of the WIMP mass. We can see that they lead, for a 1 TeV WIMP, 
to a suppression of the cross section by a factor of approximately 15. 

We now examine the differential rate with respect to recoil energy for several materials. 
This distribution is important in order to correctly apply experimental efficiencies as well 
as to assess signal-to-background levels as a function of the cut on the recoil energy, E'^*"'. 
In Fig. 1^ we present the predictions for rates differential in recoil energy on three different 
targets: Nal, ^^Ge, and ^^^Xe, including both spin-dependent and scalar contributions 
along with appropriate nuclear effects. For scalar contributions, this is the form factor 
given in Eq. (p!0|). For spin-dependent form factors, we use those presented in for 
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Figure 7: Spin and scalar nuclear form factors for Iodine as a function of the WIMP 
mass. Note that form factors are usually expressed as a function of the recoil energy 
Er = q^/2M. To generate this plot, we chose a typical recoil momentum q = /if so that 
Er{m) = m^v^M/{2{m + Mf). 



^^Ge. For Iodine, Sodium and ^'^^Xe, we adopt the spin form factors of p3[ estimated by 
considering the Nijmegen (II) nucleon-nucleon potential. 

In order to obtain the observable rates at detectors, we integrate the differential cross 
sections over the recoil energies to which they are sensitive. The minimum observable 
energy is a function of the experimental set-up and background levels. In Table ^ we 
present some current and near-future dark matter search experiments, including their 
primary target nucleus, target mass, and an estimation of the minimum recoil energy 
required for an observable event. For the Nal detectors, we have included DAMA (100 kg 
of Nal) and LIBRA (an upgrade of DAMA: 250 kg of Nal) fS^ . There are also several 
different detectors based on various isotopes of Germanium. The first stage of GENIUS 
is composed of 100 kg of ^^Ge, whereas the second stage consists of between 100-10000 kg 
of a mixture of 86% ^^Ge and 14% ^^Ge p5[. The MAJORANA experiment will search 



for double-beta-decay with 500 kg of the same mixture of 86% ^^Ge and 14% '^^Ge 
Finally, the proposed XENON experiment will consist of 1000 kg of ^^^Xe p7 . 



In Fig. P we show the potential number of events per year at detectors based on 
Germanium isotopes, assuming E^^^ is 11 keV. The bands represent potential signal 
rates as a function of the LKP mass, varying 5% < A < 15% and 115 GeV < m/j < 200 
GeV. To illustrate the importance of large mass detectors, we have chosen to compare 
the 100 kg ^^Ge GENIUS, 10^ kg ^^Ge and ^^Ge GENIUS stage II, and 500 kg ^^Ge and 
''^Ge MAJORANA experiments. Evident from the figure, in order to have even one event 
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Figure 8: Energy Spectrum of events for three types of detectors: ^^Ge (dotted curve), Nal 
(dashed curve) and ^^^Xe (sohd curve) for 171^(1) = 1 TeV, m/i=120 GeV and A = 15%. 

per year at GENIUS or MAJORANA, rrih and A must be on the small side of the band, 
and/or m^(i) must be less than about 1 TeV. In order to have ten or more events per year, 
at MAJORANA, we must have m^{i) < 700 GeV as well as small m/j and A. Thanks 
to its enormous mass, the upgraded GENIUS experiment with 10^ kg of ''^Ge and ^^Ge 
can do much better, with more than ten events per year when m^w < 700 GeV even 
for unfavorable rrih and A, and more than ten events over the entire range of m^(i) for 
optimal rrih and A. 

In Fig. |10| we show events per year, varying the parameters as above at DAMA and 
DAMA/LIBRA. DAMA, with 100 kg of Nal can observe more than ten events per year 
if the LKP is light and parameters favorable. DAMA/LIBRA, with 250 kg of Nal can 
observe the lightest relevant LKP masses even when A and rrih are at the larger end 
we consider. Finally, the XENON experiment combines a heavy ^^^Xe target with a 
large 1000 kg detector. The estimated events per year are plotted in Fig. |1T], and are 
comparable to the end-stage of GENIUS II with 10^ kg of Germanium. For comparable 
masses, XENON could in fact do better, thanks to the heavier target nucleus. 

Our results indicate that to directly detect Kaluza-Klein dark matter, heavy target 
nuclei and large mass detectors are essential. Of course, the actual reach of the experi- 
ments will depend on experimental issues such as efficiencies and backgrounds, which are 
beyond the scope of this work. However, given the relatively large event rates which are 
possible at planned experiments, further detailed study of this subject is warranted. 
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bxperiment 


iarget 


Mass 


K 


DAMA 


Nal 


100 kg 


20 keV 


DAMA/LIBRA 


Nal 


250 kg 


20 keV 


GENIUS 


73Ge 


100 kg 


11 keV 


GENIUS II 


™Ge, ^^Ge 


100-10000 kg 


11 keV 


MAJORANA 


^^Ge, ^^Ge 


500 kg 


11 keV 


XENON 


i3iXe 


1000 kg 


4 keV 



Table 2: Present and near-future dark matter detection experiments. 

6 Conclusion 

The identity of the dark matter is one of the most intriguing puzzles in modern physics, 
and has sparked a major experimental program to search for these elusive objects. In 
fact, one of the primary motivations of the supersymmetric standard model (with im- 
parity) is the fact that it has a natural WIMP candidate. However, we have recently seen 
that a large class of extra-dimensional models with a KK parity also provides a natural 
WIMP - the heavy Kaluza-Klein modes of the ordinary photon and Z boson. While KK 
parity is not a necessary feature of models with extra dimensions (just as -R-parity is an 
unnecessary feature of models with supersymmetry) , it can be imposed self-consistently in 
the low energy effective theory. Estimates of the relic density indicate that such particles 
should have masses at the TeV scale, at the frontier of current collider and direct detection 
searches. 

In this article we have made a detailed study of the direct detection of LKPs which 
scatter off of heavy nuclei. Subject to our assumptions that boundary terms are small 
(perhaps generated radiatively) and common for all quarks, our predictions depend on 
three parameters: the LKP mass, the splitting between the LKP and the first level KK 
quarks, and the mass of the zero-mode Higgs. The mixture of KK B and W3 in the LKP is 
another parameter, which would also be interesting to study in more detail. This situation 
is to be contrasted with the minimal supersymmetric standard model, which requires four 
parameters to describe the neutralino alone. We find that nuclear effects are important 
to include in the cross sections, particularly because of the heavy WIMP masses favored 
by the estimates of the relic density. These effects render rates at current dark matter 
searches small, and indicate that future experiments composed of large numbers of heavy 
nuclei can study most, but not all of the parameter space relevant for the correct WIMP 
relic density. 

Another search strategy which we have not considered here relies on indirect detection, 
in which WIMPs annihilate one another producing energetic photons and positrons [§,|11 



or neutrinos |P,|T^,|TTI] which can be observed on the Earth. The resulting LKP masses 
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Figure 9: Number of events per year for the 100 kg ^^Ge Genius experiment, the 10^ kg 
^^Ge/^^Ge GENIUS stage II, and 500 kg ^'^Ge/^^Ge MAJORANA experiment. The bands 
are obtained by varying 5% < A < 15% and 115 GeV < rrih < 200 GeV. 

which can be probed by next-generation experiments are similar to those determined here 
from planned direct searches such as GENIUS. 

In conclusion, Kaluza-Klein dark matter is well-motivated in a large class of theories 
with compact extra dimensions, and provides an interesting alternative to the standard 
neutralino LSP in a supersymmetric model. Accurate determination of scattering cross 
sections with heavy nuclei involve nuclear form factors, and indicate that future exper- 
iments can study a significant region of parameter space. At the same time, collider 
searches at the Tevatron Run II and LHC will study a similar range of parameter space. 
The exciting scenario of KK dark matter can be studied on two fronts simultaneously in 
the near future. 
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E^>20keV = 
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Figure 10: Number of events per year for the 100 kg Nal DAMA experiment, and 250 kg 
Nal DAMA/LIBRA experiment. The bands are obtained by varying 5% < A < 15% and 
115 GeV < rrifi < 200 GeV. Note that a signal in a DAMA-hke experiment is not given 
by the number of WIMP-nucleus scattering events; rather, it is the annual modulation of 
this event rate, which amounts to a small fraction of the signal rate. Experiments of this 
kind use a different methodology as they do not attempt to distinguish between signal 
and background on an event-by-event basis. Our figure should then be interpreted as a 
theorist prediction of the number of events 100 kg and 250 kg Nal detectors can detect 
per year. In the particular case of the DAMA/LIBRA experiment, the amplitude of the 
modulation must scale like the total rate. 
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Figure 11: Number of events per year for the 1000 kg ^^^Xe XENON experiment. The 
bands are obtained by varying 5% < A < 15% and 115 GeV < < 200 GeV. 
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